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Abstract 

In a number of natural and social systems, the response to an exogenous shock relaxes back to 
the average level according to a long-memory kernel ~ \/t^^^ with < ^ < 1. In the presence 
of an epidemic-like process of triggered shocks developing in a cascade of generations at or close 
to criticality, this "bare" kernel is renormalized into an even slower decaying response function 
~ \/t^^^ . Surprisingly, this means that the shorter the memory of the bare kernel (the larger 1 + 6), 
the longer the memory of the response function (the smaller 1 — 9). Here, we present a detailed 
investigation of this paradoxical behavior based on a generation-by-generation decomposition of the 
total response function, the use of Laplace transforms and of "anomalous" scaling arguments. The 
paradox is explained by the fact that the number of triggered generations grows anomalously with 
time at ~ so that the contributions of active generations up to time t more than compensate the 
shorter memory associated with a larger exponent 6. This anomalous scaling results fundamentally 
from the property that the expected waiting time is infinite for < < 1. The techniques developed 
here are also applied to the case > 1 and we find in this case that the total renormalized response 
is a constant for t < 1/(1 — n) followed by a cross-over to ~ 1/t^^^ for 1/(1 — n). 
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I. INTRODUCTION 



Many systems in the natural and social worlds are characterized by activities whose level 
A{t) at some time t is a function of its past levels {A{t), for < r < t}. This can be 
described by the generic integral equation 

A{t) = f{t)+n f A{T)^t - T)dT , (1) 

JO 

where f{t) is some source or perturbation term whose impact is instantaneous. The second 
integral term describes the propagation of past activity levels A{t) to the present time t 
mediated by the kernel ^{t — r). The summation describes that all past activities have 
an impact in the present activity level, but with a weaker and weaker weight $(t — r) as 
they recede more in the past. The "bare" kernel function must satisfy the normalization 
condition 

POO 

/ = 1 . (2) 

Jo 

Finally, the parameter n, which is in (0, 1) to ensure the absence of explosive solutions, 
describes the relative strength of triggering of future activity by past activity, as will become 



clear in the sequel. As equation ([ 
of epidemic branching models 



Ican be obtained as the statistical average of a large class 
, it is natural to refer to n as the "branching ratio" . 
We are interested in the class of systems for which the kernel expresses the existence 
of a long memory and, for the sake of concreteness, our calculations will use the specific 
form 

m - . (3) 

corresponding to $(t) ~ at long times. 

Expression ([1]) with ([3]) corresponds to a mean-field or statistical averaged description 
of the dynamics of many systems j^, such as the following examples. Present seismicity is 
in large part triggered by past seismicity over long time scales described by the Omori law 



Bfll 



% |6j, A{t) being the seismic rate in a given region above some magnitude threshold. 
Commercial and social successes have been shown to promote success over very long time 
scales . The activity A(t) here corresponds to the number of products sold or the 

number of downloads, views, attendence and so on, per unit time. Past financial volatility 
has a very long influence on future volatilty, leading to bursty intermittent behaviors also 
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characterized by long- memory power law decaying kernels 11, Q, [l^. Here the activity 
A{t) is simply a measure of the financial volatility. 
The solution of equation ([T]) has the form 

H 
10 



A{t) = f{t)+nfR{t-T)fiT)dT, (4) 
Jo 



where R{t) is the resolvent, also known at the renormalized kernel or response function, 
satisfying the equation 

R(t) = <^{t) +n f' <!>{t-T)R{r)dT . (5) 
Jo 

The name "response function" refers to the fact that, if f{t) = M6{t — to) is a delta function 
describing a sudden impulse of amplitude M applied to the system at some time to then, 

A{t) = M nR{t- to) (6) 

for any t > to. The resolvent R{t) thus describes the response of the system to an impulsive 
perturbation. 

The focus of the present paper is to analyze the properties of the response function R{t) 
by addressing two specific questions. 

• PARADOX 1: at criticality, the shorter the memory of the bare kernel, the longer the 
memory of the response function! 



For the c 



ass of systems with long memory with < ^ < 1, starting with MontroU and 



Scher a number of authors have shown that 

R{t) ~ 1/t^-^ forn = 1 or if n < 1, for t < t* , (7) 

where t* = p [^^-^^f^^ . For t > t*, R{t) ~ <l>(t) ~ l/t^+^ Q, H, Q (see in 
particular 17| for a synthesis). The value n = 1 corresponds to the critical regime. 
While the derivation of this result is rather straightforward when using the Laplace 
transform operator as we recall below, this result is paradoxical. Indeed, at face value, 
it means that the larger 6 is, the shorter is the memory encoded by the bare kernel 
and the longer is the memory described by the response function R{t). A first 
goal of this paper is to solve this paradox by a detailed analysis of the role played by 
the cascade of triggering events intrinsically embodied in equation ([1]) for n = 1. 
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PROBLEM 2: power law exponent of the resolvent for 1 < 6 < 2 at criticality (n = 1). 

Recently, an empirical example of the regime where 1 < 6 < 2 has been discovered 
in the humanitarian response to the destruction brought by the tsunami generated by 



the Sumatra earthquake of December 26, 2004, as measured by donations [18|. The 
data suggests that n < 1 so that the observed response is R{t) ~ ~ l/t^~^^, i.e., 
the effect of multiple triggering occurring for n ~ 1 did not occurred in this episode. 
Here, we ask what would be the response function of a such a system in which the 
bare kernel is of the form ([3]) with 1 < < 2 if an epidemic cascade characterized by 
n ~ 1 occurred. Clearly, the solution ([7]) cannot apply for 1 < 6 < 2, since it would 
lead to a growth of the response function with time! The solution for this regime that 
we provide below further illuminates the solution of Paradox 1. 

In a nutshell, the solution of the paradox developed in the sequel of this paper is based 
on the decomposition of the total activity as the sum over a time- varying number K{t) ~ t^ 
of generations which have been activited until time t. Since each generation k, with 1 < A; < 
K(t), contributes to the total activity with an amplitude which is proportional to ~ kn^ /t^^^ , 
the total activity is therefore ~ jitf x (^fi? kn^) ~ [K{t)f /t'^+'^ ~ t^e/^i+e _ 
Therefore, the larger 6 is, the shorter the memory of but the faster growmg IS 

the number of generations that are triggered up to time t. Thus, the resulting slower 
decaying renormalized response function ~ results from the fact that the number 

of triggered generations grows sufficiently fast so as to more than compensate the shorter 
memory associated with a larger 9. This anomalous scaling results fundamentally from 
the property that the expected waiting time is infinite for < < 1. Adaptation of 
this reasoning to the case 1 < 9 implies that the renormalized activity is a constant for 
t < 1/(1 — n) followed by a cross-over to ~ Ijt^^^ for t ^ 1/(1 — n). 

The rest of the paper develops the derivations of these results and is organized in the 
three following sections. The next section 2 constructs the mathematical building blocks 
used in the subsequent sections. In particular, subsection III CI summarizes the main results 
obtained using the Laplace transform applied to expression (j5]). Section 3 presents a detailed 
derivation of the solution of Paradox 1 in terms of a generation-by-generation decomposition. 
Specifically, subsection IIIIBI dissects the two key contributions in the global activity and 
subsection IIII CI presents the intuitive derivation of Paradox 1 based on scaling arguments 
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for long waiting times. Section 4 applies the same approach in terms of the generation-by- 
generation decomposition for the case 1 < 6 < 2. 



II. MATHEMATICAL BUILDING BLOCKS 

A. General relations obtained by using the Laplace transform operator 

The standard way to solve the integral equation ([5]) is to apply the Laplace transform. 
This transforms the integral equation into an algebraic one 

R{s) = ^{s) + n^{s)R{s) , (8) 

where the tilde denotes that the corresponding function is the Laplace image of the original 
function. For instance, the Laplace image of the kernel is 

POO 

$(s) = / (^{t)e-''dt . (9) 
Jo 

It follows from ([8]) that the Laplace image of the resolvent is given by 
Accordingly, one can represent the solution A{t) of ([6]) in the form 

oo 

A{t) = Y.A,{t) , (11) 

k=l 

where 

Auit) = Mn^^uit) (12) 

and $fc(t) is the inverse Laplace image of (^^{s). 

The series (111 I) has a transparent meaning, when interpreted in the context of epidemic 
processes. Interpreting f{t) defined in (j4]) as some ancestor event of amplitude M, then A{t) 
is the mean birth rate of its offsprings. Correspondingly, the fc-th term Ak{t) in the series 
(|TT1) is the mean birth rate of those offsprings of the fc-th generation. In this context, n is the 
critical parameter of the corresponding branching process. In what follows, it is convenient to 
use the terminology of branching processes in order to describe the characteristic properties 
of the resolvent R{t) and the corresponding solution A{t) of the integral equation ([1]). For 
simplicity, but without loss of generality, we put M = 1. 



B. Scaled kernel and resolvent 



Specific calculations will be performed with the form ([3]) of the kernel. But, whenever 
possible, we will keep the discussion as general as possible. In particular, we will consider 
the general class of kernels of the integral equation ([T]) which has the form 

4W = i^(i), (13) 

where ^ > is some unique characteristic time scale, while the kernel tf{x) is some given 
function of its dimensionless argument, satisfying to normalization condition 

roo 

/ !f{x)dx = 1 . (14) 
Jo 

For the choice ([3]), we have 

^(^) = (^^T)m • (15) 
We will restrict our study to the case of kernels possessing the power asymptotics ~ 
t~^~^ for t oo, with < 6 < 2. The case 6 = 1 requires a special treatment. We will 
not present it for the sake of conciseness, while it is clear that the method presented below 
allows one to easily provide the needed detailed description in this case. We note that the 
main scaling laws obtained below remain valid for 6 = 1, while some details and corrections 
to scaling differ. 

For the family of kernels given by f|T3|) with a single characteristic scale g, the resolvent 
R{t) and the mean activity A(t) can be represented in forms analogous to (IT^ : 

R(t) = -nl-,n] , A(t) = -Al-,n] , A{x , n) = nn{x , n) . (16) 
Q \g J Q \Q J 

The Laplace images of the resolvent TZ{x, n) and of the mean activity A{x, n) are given by 
ny, n) = ^^y] = - f; n^^\y) , A{y, n) = n'^\y) , (17) 

where 

POO 

(p{y) = / ip{x)e-''ydx (18) 

JO 

is the Laplace image of the bare kernel '{>{x). In particular, the Laplace image of the kernel 
given by (ITSll is equal to 

$(s) = 6 {qsY e^' V{-6, Qs) , ^ (p{y) = 6y'eyV{-6, y) . (19) 
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The asymptotic power law of the kernel ip{x) ^ x ^ ^ (x ^ oo) leads to the following 
asymptotics for the Laplace image (p{y) for small y values: 

(p^y) c^l + ay-Py' , < 1 • (20) 

For the particular case (IT^ . we obtain 

^(y)^i + ^-y^r(i-0) , (21) 

In this case 

a ^ a{d) = , m = r(l - e) . (22) 

Notice that a{9) and (3{9) change sign as 9 crosses the value 9 = 1. Specifically, a{9) and 
f3{9) are negative for 9 G (1, 2) and positive for 9 G (0, 1). These signs are the consequence 
of the fact that $(t) given by (fT5|) is one-sided, i.e., identically equal to zero for any t < 0. In 
the following, we will thus consider the general asymptotic expression (|20|) with coefficients 
a and /? of the same sign. 

C. Direct derivation of the resolvent R{t) 

Substituting fl20|) into (|T0|) and using its first equality leads to 

^(^'^)- ^ q = l-n. (23) 

q-ay + (3y^ 

1. Case < 6* < 1 

In this case, the term ay can be neglected in (|23|1 . which becomes 



We have intentionally replaced (3 by \(3\ to stress the important fact that, for 9 G (0, 1), the 
parameter (3 is positive as can be checked from its explicit value given in ( l22l) . For very 
small y, expression (IMl) can be further expanded into the following asymptotic relation 

ny,n)^^--^y\ (25) 
q q 

Using the standard correspondence between functions and their Laplace transforms 



we obtain the asymptotic time dependence of the resolvent: 

7^(x,r^)^^x-^-^ x » ^ . (27) 

The intermediate asymptotic regime corresponds to an interval of still small values of but 
not too small so that the following inequality \l3\y^ ^ q holds. Then, one can neglect the term 
q in the denominator of expression fl2^ to obtain the following intermediate asymptotics 



^(^'^)^^ - 7^(x,r^)^^-^a:-^+^ (28) 

The two regimes (!27|) and (l28ll are asymptotics of the inverse Laplace transform of (12^ 
whose explicit expression reads [l^ 

where 

sin(7r6') n , /""^ u^e~^du ,„ , 



For instance, for = 1/2, 



Q(x, 1/2) = ^ - e"erfc(v^) . (31) 



2. Case l<9 <2 
We rewrite (|23|) as 

For 



^(y,^)^ ^1 / , • (32) 
q+ \a\y- \l3\y^ 

1 



ayl < 1 , (33) 



then 



^(y,n)^i-M^ + M/, (34) 



q q^ q^ 



which is the Laplace image of 



7^(x, n) ^ -M- X-'-' = 4 X-'-' , X » i . (35) 
The intermediate asymptotic describes the time interval such that q <^ \a\y, leading to fl5^ 

'^(^•"'"|ato-W = J-|a|-L-'- 
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For sufficiently small \y\, one gets 

n{y,n) 



fI y 



+ 



/3 



■ y 



(37) 



which is the Laplace image of 



7l{x, n) 



1 



/5 



X 



(38) 



T{2-e) q 

Expressions ( !38l) and (|35l) show that the resolvent TZ{x, n) is a constant plus a weak 
power law correction ~ l/x^~^ for 1 <^ x ^ 1/(1 — n) which crosses over to ~ l/x~^~^ for 
X ^ 1/(1 — n). A function proportional to the resolvent is plotted in Fig. 4 at the end of 
the paper, which shows these different regimes (138|1 and (135|1 . 



D. Asymptotic of the mean activity 

While the previous subsection provides the expressions of the different regimes of the 
resolvent, their derivation using the Laplace transform augmented by different expansions 
do not provide an understanding of the derived terms, which we would wish to be based on 
the underlying mechanism of cascades of triggering over different generations. To achieve 
this goal and remove Paradox 1, we have to explore the asymptotic behavior of the mean 
activity A{x, n) for large x values. For this, we study the asymptotics of the corresponding 
Laplace image Aiy, n) for small y values. Substituting in the series (fT7|l the asymptotic 
expression fl20|) and using the asymptotic relation 

{I + ay - [5y'^f ^ e^'^y-^^y', |a|/-/?/|<l, (39) 

we obtain 

oo 

A{y, n) ~ 5] n^i){\(3\^'^k^'%- ^)e°'=^ , (40) 
fc=i 

where 

^{y-e) = e-'''^'^^f^^y\ (41) 

Taking the inverse Laplace transform of the series fHOj) . we obtain a series representation of 
the sought mean activity as 

CO 

A{x,n)^Y.^k{x,n) , (42) 

k=l 
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where 

M^. n) = ^17^^ (f7^' ^) ' (43) 
and ipi^] ^) is a stable distribution, whose Laplace image is given by expression ( 14T|) . 

The asymptotic validity of relation (H2l) for the dependence of ^(x, n) for a; ^ 1, describ- 
ing in particular the case where the kernel is given by expression (fT5l) . transforms into an 
exact equality for the mean activity if the kernel of the integral equation ([1]) coincides with 
the stable distribution 

^t) = -^p(-■,e] , (44) 
whose Laplace image is given by fHTl) . 



E. Properties of the stable distribution 'ip{x;9) defined by (|41|) and (I44p 

Formulas 0421) and fl43p imply that a better understanding of the asymptotic shape of the 
mean activity A{x, n) is dependent on a detailed knowledge of the properties of the stable 
distribution il){x]9). This subsection is devoted to this question. 

There are many integral representations of the stable distribution 'il){x\ 9). In particular, 
one can show that 



ij{x-e) 



TT Jo 



exp 



cos 



U X 



(45) 



cos \ux + u sm I — I sign(^ 



1) du 



o<e<2 , e^i . 



Explicit analytic expressions of the stable distribution if){x; 9) exist for some specific values 
of the parameter 0. For illustrative purposes, we will use below two such stable distributions. 
The first one is the famous Levy stable law 



ij{x- 1/2) = exp (-^ ) , {6 = 1/2) 

IX\/tix \ Ax, 



(46) 



and the other is 



^(x;3/2) 



TT 



5 2 4x^ 



6' 3' 27 



'7 4 4x^' 



6' 3 27 



(47) 



{9 = 3/2) 



where iFi (a, b, c) denote a confluent hypergeometric function of the first kind. 
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All stable laws possess in common long and short tails. By definition, the long tail is 
their power law behavior at a; ^ oo: 

il){x] 9) ~ , X oo . (48) 



These short tail of stable distributions consists in a very fast decay of ip{x; 9) to zero as 
x— >-OforO<^<l corresponding to an essential singularity at x = 0, and in their super- 
exponentially fast decay as x — oo for 1 < 9 < 2. The following asymptotic formula is 
true 

2-e / , , e - 

9-1 



^(x; 9) ^ ^ ^ ( M ) exp { -\9 - V ' 



(49) 

x^0+, if O<0<1, 

X -oo , if I <9 <2 . 

Curiously, for 9 = 1/2, this asymptotic formula coincides with the Levy stable law 
According to the asymptotic relation fH9l) . the stable distribution ■ip{x]3/2) decays at x 
— oo according to 



^(x;3/2) ~ y^exp (^-^|x|^^ , x ^ -oo . (50) 

In practice, the asymptotic expression such as flSUl) can be verified to be extremely accurate 
already for x < — 2. 

III. SOLUTION OF PARADOX 1 FOR < 6* < 1 

The asymptotic behavior for x ^ 1 of the mean activity A{x, n) given by fH2]) is qualita- 
tively different for 9 G (0, 1) and for 9 G (1, 2). In this section, we focus on the former case 
^^G(0,1). 

A. Integral approximation of the mean activity A{x^ n) 

As summarized in the statement of Paradox 1, for 9 G (0, 1) and close to criticality 
(n < 1), the mean activity A{x, n) exhibits a double power law behavior, with the coexistence 
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of the power law asymptotic 

A{x,n)^x-^'\ x>l. (51) 

for very large x values, and of an intermediate asymptotic regime for smaller x (but still 
remaining large) 

A{x,n) ~ x^^^^ , intermediate asympotics . (52) 

The goal of this subsection is to show that the intermediate power asymptotic fl52|) results 
from the fast decay of the short tail part of the stable distribution iIj{x;6), which controls 
the mean activity Ak{x, n) defined by ( l43l) of the activity resulting from the k-i\i generation. 
In contrast, we will show that the power asymptotic (15T1) is due to the corresponding power 
asymptotic of the bare kernel $(t) ~ t~^~^. 

The first step consists in noting that, for < 6 < 1, the shift in expression fH3l) for 
Ak{x,n) can be written as 

Oik , 0-1 , 

k— . (53) 



\f3\i/ef.i/9 

It thus tends to zero for A; oo, so that one may neglect it as it will not impact the 
asymptotic law of the mean activity A{x,n) for large x's. In other words, for < < 1, 
one may without essential error replace the k-th generation mean activity A{x, n) by 

fx \ 

M^, n) = ^TJe^e^ [w^F'^ ^) " ^^^^ 

The next step is to notice that, if x large enough, then Ak{x,n) becomes a sufficiently 
smooth function of the argument k and one may, without essential error, replace the series 
fj42|) by the integral 

A{x, n)c^J Ak{x, n)dk ~ I j^TJe]^'^ ^) ' ^^^^ 

in which we have defined 

, = . (56) 

We suppose everywhere below that I7I ^ 1 (n is close to its critical value 1), so that the 
exponential function e~"'^ in the integral (l55l) is also a smooth function of k. 
Using the change of variable 

_ X L _ ^ (^\^ ^ fxY du 

""'W^ ^ MW ' mUJ 17' ^^^^ 
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the integral (l55l) becomes 



where 



As ^ 1 for X ^ 1, one may replace without essential error the upper limit in the 
integral ( l58l) by infinity: 

For 7 = = 1), we obtain the mean activity given by fl60l) as the power law ( 152|) 
corresponding to the intermediate asymptotics: 

Ant(x,n)^^p, (61) 

where the index 'int' refers to "intermediate asymptotics" and 

(62) 

It is reasonable to choose (3 = r(l — 6) defined in (122|) and using the well-known identity 
r(l - 9)T{1 + 9) = 7T9csc{n9), we obtain 

, , , sin(7r^^) -I ,n 

TT 

For 7 = (ra = 1), the intermediate asymptotic fl63|) invades the whole large x regime, so 
that its determination is clearer. 

The fact that expression ( l63l) holds only for an intermediate range of x's values for 7 > 
{n < 1) is now determined from the following derivation. Let C,{9) be such that, for if 
X > C,{9), the stable distribution ip{x\9) is not different from its long tail (l48l) within a 
specified error margin. For example, for the Levy stable distribution il){x\l/2) given by 
l6|) . ^{9) can be taken equal to 3 when considering an error margin of less than 1%. 
From expression ( l60i) . one can see that, if the following condition holds. 



7 / X 



>i ^ ^^m[ — ] (64) 



1^1 \m) ~ ' V7 
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then, without essential error, the stable distribution ip{x; 6) can be replaced by its long tail 
in the integral fl60l) . This leads to the following approximate relation 

, 1 9 [°° ( 1 du 



Changing the integration variable to 2; = u/x, we obtain 

1 OP {6) 



^(x,n)^--— (66) 



where 

B. Dissecting the two key contributions to solve Paradox 1 

The asymptotic analysis of the mean activity A{x, n) presented in the previous sub- 
section, while sufficiently rigorous, does not provide an intuitive understanding of the two 
intermediate and asymptotic regimes and of the solution of Paradox 1. In the present subsec- 
tion, we provide a cruder but more transparent analysis, which reveals the hidden springs 
of the crossover from the intermediate asymptotic power law (!52l) to the full asymptotic 
power law( l5T]) . The next section will provide a different approach which illuminates even 
further the mechanism of the transition from the bare kernel time decay to the resolvent 
time dependence. 

For this, we replace in expression (j^3l) the stable distribution ipi^x; 6) by its "geometrical 
power law skeleton" 

M^;e) = \ (68) 



, x< m . 

Replacing in fll51) ip{x]6) by ipQ^x; 6) (and neglecting the shift parameter flS5]) ). we obtain 

Al{x, n) ^ - W"k^"m) , (69) 

where 1(2;) is the unit step function. This term l(x — \l3\^^^k'^^^^{6)) is important, as it 
accounts semi-quantitatively for the fast decaying short tail of the stable distribution ip{x] 9), 
which has the role of effectively truncating the series fl42p at large fc's. 
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For concreteness, let us consider the particular kernel defined by (fTSi) . for which 
(3 = r(l — 6). We can then rewrite relation fl69|) in the more transparent form 

Al{x, n)=e x-^-^ k l{k{x, e)-k) , (70) 

where 

x^ x^ 

" " e'Wr(i-e) ■ ^^^^ 

Choosing for simplicity ^^{0) = we obtain 

k{x, e)=x^ . (72) 

Substituting expression flTOj) into the series P3|) . we obtain the mean activity estimated in 
this geometrical skeleton approximation, denoted as v4°(x, n), as 

^°(x, n)^e x-i-^ S{k{x, e),n) = e ■ 5(x^ n) , (73) 

where 

S{K,n) = J2k-n^ . (74) 

fc=i 

Expression (1731) with (I74p allows us to pinpoint the origin of the slow decay ~ of 
the resolvent in the intermediate asymptotic regime or for n = 1 as due to the fight between 
the fast decay ~ of the bare kernel and the growth ~ of the contributions to 

the activity at time t of all generations set in motion up to time t. This later growth is 
controlled by the short tail of the stable distribution iIj{x; 9) corresponding to the truncation 
on the above geometrical skeleton approximation fl68l) . The main contributions to A{x,n) 
are provided by the first k{x,6) summands, because the mean activities Ak{x,n) of the 
highest order generations, for which k > k{x, 6), are, for given x = t/ g, not yet large enough 
to influence significantly the total activity level A{t). Roughly speaking, the larger the order 
/c of a generation, the later its contribution Ak{x, n) is felt. 

For n = 1, the sum (1741) reduces to 



Using (1721) . we obtain 



S{K,l) = j:k = \K{K+l)^\K' . (75) 



s{k{x,e)^)^\x^' . (76) 
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Substituting this expression (1761) into (1731) . we finally obtain promised power law (152|) 

A\x, n) ^ ^ ■ X-'-' -x^'^^-- x-'+' . (77) 

This last equation (1771) illuminates the origin of Paradox 1: the larger 6 is, the faster the 
decay of the bare kernel ~ but the larger the number k{t, 6) ~ of generations which 

are activated up to time t and the greater their combined contribution ~ [k{t, 6)]"^ ~ t"^^ to 
the overall activity at time t, so that, all being taken into account, the response function 
actually develops a longer memory ~ x = Paradox 1 can thus be seen as 

a result of an "anomalously" slow triggering of successive generations associated with the 
infinite average waiting time between triggered events. Indeed, the average waiting time 
between two events, as described by the bare kernel, is defined by 

(t) ~ hmr^+oo J -^dt , (78) 

which is diverging as the upper bound T of the integral goes to infinity, for < 1. This 
divergence is a standard diagnostic of the existence of an anomalous trapping time regime 
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2l| . leading to anomalous scaling laws. In the present case, the "anomalous" scaling 
law is the "renormalization" of the bare kernel time decay ~ into the resolvent time 

decay ~ The next subsection IIII CI re-derives this result from scratch by using a 

completely intuitive and straightforward reasoning, exemplifying that the root of Paradox 1 
indeed lies on the diverging mean waiting time between triggered events and its associated 
anomalous diffusion. 

But before doing so, we exploit the present analysis to describe the subcritical case n < 1. 
For arbitrary n, the sum ( l74l) is equal to 

S{K,n) = ^[1 - {1 - qni + K - q)] , q=l-n. (79) 



For small enough g, such that Kg <C 1, expression (179|) has the same square asymptotic (175|) 
as for n = 1. In contrast, for nq ^ 1, S{K,n) tends to a constant limit S{K,n) — > ^. This 
implies that, when k{x,9)q ^ 1, the mean activity A^{x,n) obeys the same power law as 
the bare kernel, 

A\x,n)^^-x-'-\ (80) 

as seen from (173|) . The function S{K,n) defined by (179|) is plotted in Fig. 1 for different 
g = 1 — ?7. values. It demonstrates the crossover from S{K,,n) ~ for Kg ^ 1 to S{K,n) ~ 
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const for Kq ^ 1. This crossover governs the crossover of the resolvent from ^ for 
t < 1/(1 - n)!/^ to l/t'+^ for t > 1/(1 - ny/^. 
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CO 
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I 

Fig. 1: Plot of function S{K,n) defined by ([79]) for g = 1 - n = 0.01; 0.015; 0.02, 
illustrating the cross-over between the power quadratic behavior (j75p for Kg ^ 1 to a 
constant for Kq ^ 1. 

C. Intuitive derivation explaining Paradox 1 based on scaling arguments for long 
waiting times 

Let us now re-derive all the key results of the previous sub-section by a completely 
different and intuitive route. Our approach is based on the conceptual view of the total 
activity A{t) at a given time t as the superposition of the activities Ak{t) coming from all 
possible generations k = 1,2, ... that are significant at this time t. 

Consider a first burst of activity starting at time 0, constituting the event of zero-th 
generation. This initial event may lead to an event of first generation at a latter time ti, 
which itself may trigger an event of second generation at time t2, and so on. We assume 
that these events constitute the starting time for each successive generation to contribute 
significantly to the overall activity. For a given time t of observation of the activity, our 
problem is to determine the typical time tk{t) of occurrence of the k-th generation and its 
corresponding contribution Ak(t). 
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K{t) X / / dt/t^+' ~ 1 . (82) 



For the first question, we use the interpretation that the bare kernel is nothing but 
the probability density function (pdf) of the waiting time t from a burst of activity and its 
first triggered event. Let us suppose that K(t) generations have been triggered over the 
total time interval t. Time consistency imposes that 

ti+t2 + ...+tK(t) =t . (81) 

Let us call tmax(i) the largest waiting time among the K(t) values ti,t2, ...,tK{t)- Since the 
probability that a waiting time between two successive generations is equal to or larger than 
^max(i) is of the order of Itl^{t) dt/t^^^, by consistency, one must have 

n+OO 
' ''max 

it) 

Expression (15^ just states that there is typically just one waiting time of the order of 
the maximum waiting time tmax(^) among the K{t) waiting times between the successive 
generations. The solution of (152]) is 

tmax(t) ~ Q[K{t)]'/' . (83) 

The dependence K{t) as a function of t is then obtained by estimating the l.h.s of (!8T|) as 

„ / ''max 

ti + t2 + ... + tKit) = Kit) X {t)t ~ Kit) X dt t/t^+' ~ Q[Kit)f' . (84) 

We have used the fact that the average waiting time {t)t between successive generations 
has to be estimated by the standard sum of all possible t's weighted by their corresponding 
probability, but with an upper bound since no waiting times larger than tmax(i^) are sampled 
in the finite set of K{t) realizations. This trick is standard to tame the infinities of the 
unconditional average waiting time (t) defined by (!78|) in the limit T — >■ +oo, found for 
6 < 1, leading to anomalous diffusion and other abnormal scaling effects 
(IHTj) . we obtain [K{t)Y/^ ~ t/g and thus 
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2l|. Then, by 



K{t) ~ (t/g)' , (85) 

which retrieves f l72]) obtained in the previous sub-section. 

The contribution Ak{t) of the k-th active generation at time t has two important terms. 
The first one is the probability n'' that k generations have occurred. The second one is 
based on the concept that each activated generation contributes proportionally to the bare 
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kernel ~ 1/ 1^^^ but with a characteristic time scale pk equal to the only existing time scale 
associated with the generation k, namely its waiting time t^: pfc ~ tk until its happenance. 
In complete analogy with the shape of given by ([3]), this leads finally to 



(86) 

Ordering the indices by increasing values of the waiting times tfc(t), expression ( l85l) implies 
that tfc(t) ~ gk^^^, and thus 

ek 



The total activity is thus 



^fe ~ ^ ^ ' for t > tk{t) . (87) 



K(t) n K{t) 

At) = E Mt) -TTTe^Y.kn\ (88) 



fc=l k=l 



which recovers (1751) with flTll) . 



D. Validation with the exactly solvable case 6 = 1/2 



It is always useful to check the validity of asymptotic relations when exact results are 
available. Here, we compare the above asymptotic relations of subsection IIII Bl for the mean 
activity A{x,n) with the series f H2l) . fH3l) for 6' = 1/2. In this case, the series ( H2l) takes the 
form 

(89) 



, . V 1 ( X -\- 2/c 1 



k=l 



while its geometrical skeleton is equal to 



1 



n) . 



(90) 



A more accurate approximation of the series (15^ than that given by integral (jSO]) in the 
case 6 = 1/2 is equal to 



Ait(a;,ri) 



1 



1 — 7a/x exp i j erfc 



(91) 



where 7 is defined by f l56|) . For = 1/2, expressions fl6Tl) and fl66|) become 

1 



Ant(a;,n) 



1 



7 = 



(92) 



_ 2xi/x7^ ' ~ 7^ 
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Fig. 2 compares the mean activity A{x,n) given by (189|1 . its geometric skeleton ^''(x, n) 
obtained from fl90|) and its integral approximation Aint{x,n) given by fl9T]) . One can check 
that they are practically undistinguishable for any x > 10. One can also observe the two 
power law regimes A ~ x^^^^ and their cross-over. 




X 



Fig. 2: Plots of the mean activity A{x,n) given by (j89p . its geometric skeleton 
A^{x,n) obtained from (f90l) and its integral approximation Amt{x,n) given by (i9T|) . 
for = 1/2 and 7~g = 1 — n = 0.01. The dotted straight lines show the limiting 
power asymptotics (|92]) corresponding to ^ ~ x^^^^. 



IV. PROBLEM 2: POWER LAW EXPONENT FOR THE RESOLVENT FOR 1 < 9 
A. Integral representation 

It is convenient to re-express equations (l42l) and ( l43l) in a form more adapted to the case 

e>i: 

For large x ^ |q;|, one may, without significant error, replace the series in (1931) by the 
integral 

/oo 
e-^'' Ak{x-e)dk , (94) 
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where 

The "body" of the stable distribution 6) is concentrated near u = 0, which means that 
the "body" of Ak{x; 9) taken as a function of k is concentrated in the vicinity of 

K = ^. (96) 

\a\ 



As a consequence, the value of the integral ( jMl) is qualitatively different depending on the 
value of the parameter 

£ = ^fc, = ^ . (97) 
|a| 

B. Early time asymptotic 

Let us first consider the regime £ ^ 1. In this case, one may, without significant error, 
put 7 = in (|94|) to obtain the following approximate relation 

/oo 
Ak{x;e)dk . (98) 

Given the definition fl56l) . putting 7 = is equivalent to neglecting the difference between 
the branching ratio n and its critical value 1. In other words, for times sufficiently short 
such that e defined by (1971) is small, the mean activity is faithfully described as if the system 
was in the critical regime n = 1. Taking into account that the effective width of the function 
Ak{x; 9) of the argument k, defined as the domain in which Ak{x\ 9) is significantly different 
from zero, is much smaller than fc*, one can replace k by the constant k^ given by fl96l) in the 
denominators of the r.h.s. of expression (195|) . Then, using the change of integration variable 
k u = ^^!"^i/o ) we rewrite the integral (!94l) in the approximate form 

Aint{x,n) ~ — / ^(m; 9)du = — (^ < 1) , (99) 

\a\ J-00 \a\ 

as result of the normalization condition of the stable distribution tplx; 9). In sum, we have 

^(x,n)~^ = ^-l , ^«1 ^ a:« J. -, (100) 

|a| \a\ [9 — — n) 

where we have used definition fl22l) for a and fl56l) for 7, assuming that n is close to 1 so that 
7 ^ 1 — n. The result f llOOp expresses that, for n ~ 1, there is plateau in the mean activity 
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A{x, n) as a function time, for early time x <^ 



1 



. As n moves closer and closer to 



(e_i){i_„) 



1, the regime where A{x,n) ~ — 1 extends to longer and longer times. We note that this 
constancy of the resolvent at criticality n = 1 is well-known for an exponentially decaying 
kernel function ~ exp[— rt], corresponding to the resolvent R{t) ~ exp[— r(l — n)t]. 
The novel behavior found here is the existence of a non-trivial cross-over to the power law 
~ l/t^~^^, as explained in the next subsection IIV C[ 

This result (llOOp can be recovered simply by using the arguments of subsection IIII CI 
They extends straightforwardly to the case 6' > 1, for which the unconditional mean waiting 
time (t) defined by ( 1781) in the limit T — > +oo is now finite and well-behaved. This implies 
that ti + t2 + ... + tK(t) in (IHTl) is well approximated by K{t) ■ (t), which, being equal to 
t, yields K(t) ~ t, for all values of 6* > 1. This extends the result ( 185|) previously derived 
only for < < 1. Then, expression fl87|) still holds and we obtain finally that expression 
(IHHj) holds with K{t) ~ t. For n = 1, we recover that the leading term describing the time 
dependence of A(t) is a constant as described by ( llOOp for the case 9 > 1. The simple scaling 
argument of subsection IIII CI shows that this result can be actually extended to all positive 
values of 6 > 1. For n < 1 but 1 — n small, A{t) ~ constant still holds for times ^ 1/(1—72). 

We can combine the results of the intermediate asymptotics valid for x <^ {e-i){i-n) 
6 G (0, 1) and for 6 > 1 hj the following power law dependence of the mean activity 



A{x, n) ~ X 



(101) 



where the exponent 5{9) is given by 



i-e 



< ^ < 1 , 



5{e) = {I - e)i{i - 6) 



(102) 







^ > 1 . 



Fig. 3 shows the exponent 5{6) given by (11021) as a function of 6. 
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Fig. 3: Dependence of the exponent 5{9) given by (jl02p of the mean activity A{x,n) 
as a function of 0, in the intermediate power asymptotic defined 1 — n <^ 1. 



C. Long time asymptotic 

Let us now consider the asymptotic behavior of integral (|94l) for £ ^ 1, where e is defined 
in ( 1971) . In this case, the main contribution of the integral is the power law ( l48l) of the stable 
distribution, so that one may replace the function Ak{x; 6) by 



|r(_^)| {x-\a\ky+^ 
Accordingly, integral (IMl) takes the form 

rK-i k e-^'' dk 



(103) 



Ant(a;,n) 



^ > 1 • 



(104) 



\T{-e)\ Ji {x-\a\ky+^ \a\ 

The upper limit of this integral removes the influence of the singularity at k^, defined by f l96|) 
which is irrelevant for e ^ 1. One may interpret the upper limit as resulting from the short 
tail of the stable distribution. 

Due to the fast decaying exponential e"'^'^, fll04p can be approximated by 



T{-9) 



X 



r k e-^'' dk e = ^ > 1 , 
Jo \a\ 



(105) 
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which finally leads to 



Amt{x,n) 



p 



X 



-0-1 



If the kernel is of the form ([T5]), then P = T{1 



e 



^ — X 

T 



e > 1 , 
) and one has 
£ > 1 . 



(106) 



(107) 



D. Exact results for 8 = 3/2 



It is useful to check these results for 6 = 3/2, for which we can make use of the explicit 
expression ( l47l) . To be specific, we also assume that the kernel $(t) is given by formula ( fTSl) . 
so that the parameters \a\ and |/3| are equal to 



2v^ 



\a\ 



Accordingly, the series in (19311 reads 

oo 

A{x,n) = J2 



n 



k=l 



x-2k 



3/2 



(109) 



while the power asymptotics fllOOl) and (11041) take the form 

' 1 7X 
2 ' 

A{x,n) ~ 



2 

7X 



«1, 
> 1 . 



(110) 



I 272 x5/2 ' 2 

Fig. 4 plots A{x, n) given by f ll09p and its asymptotics given by flllOp as a function the 
reduced time x for n = 0.99. One can observe the predicted constant plateau for times up 
to ^ 1/(1 — n) ^ 10^, followed by the power law ~ l/t^~^^ at larger times. 
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X 



Fig. 4: Dependence of expression (|109p for the mean activity (solid line) and its two 
asymptotic regimes (jllOp (dotted lines), for 9 = 3/2. The branching ratio is equal to 
n = 0.99, i.e., 7 = 0.01. 
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